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ABSTRACT 

The environs of three supernova remnants (SNR) with associated OH 1720 
MHz masers, G349. 7+0.2, CTB 37A and G16. 7+0.1, have been surveyed in 
the CO J=l-0 transition with the 12 Meter Telescope of the NRAO, using the 
On-The-Fly technique. These observations have revealed a number of molecular 
clouds interacting with the SNR shock fronts. Most of the OH 1720 MHz masers 
have been found to lie over CO concentrations, and the maser velocities are 
coincident with the CO peak velocities to an accuracy better than 2 km/s. The 
present data trace the interstellar medium (ISM) structures interacting with 
the SNRs; however, to probe the shocked molecular gas in which the OH 1720 
MHz emission originates, higher excitation transitions and more complex species 
should be observed. In CTB 37A, where the shock velocity into the molecular 
cloud could be determined, it has been found to be of C-type, in agreement with 
theoretical predictions. Part of the rim of G16. 7+0.1 appears to be flattened by 
a dense external cloud, yet the only associated OH 1720 MHz maser lies near 
the opposite region of the remnant. This behavior, also observed in IC 443 and 
3C 391, seems to contradict the suggestion that OH 1720 MHz maser emission 
occurs mainly for transverse shocks. 

Subject headings: ISM: clouds — ISM: molecules — masers — ISM: individual 
(G349. 7+0.2, CTB 37A, G16. 7+0.1) — supernova remnants 
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1. Introduction: 

When a SNR expands near or inside a molecular cloud, the shock driven into the cloud 
can accelerate relativistic particles, heat and compress the molecular gas, and change its 
chemistry, in addition to being an important source of turbulent mixing. If the turbulence 
is not strong enough as to disrupt the cloud, condensed clumps may be created which 
eventually may become new stars. Therefore, the interaction of SNRs with molecular clouds 
joins the two endpoints of the life cycle of stars and gas in our Galaxy (Frail, Goss & Slysh 
1994). The study of such interactions sheds light on several different fields in astrophysics. 

Unfortunately, confusion introduced by unrelated gas along the line of sight, makes 
it not trivial to establish unambiguously whether an SNR is physically associated with a 
molecular cloud or whether their distances differ and they appear as positionally coincident 
just because of a projection effect. Morphological signatures like arcs of gas outlining 
an SNR, or indentations in SNR edges surrounding gas concentrations, are often used 
(Landecker et al. 1989, Reynoso et al. 1995, Dubner et al. 1999). More convincing, but 
also more rare, are line broadenings (e.g. Frail & Mitchell 1998, Reach & Rho 1999) or 
wings (Seta et al. 1998) indicative of shocked gas. 

Almost three decades ago, Goss (1968) discovered for the first time OH emission in 
the satellite line at 1720 MHz associated with two SNRs: W28 and W44 (see also Goss & 
Robinson 1968). Hardebeck (1971) found that this emission comes from several compact 
regions with strong inversion of the OH levels. Later, another OH 1720 MHz emission 
feature was found towards IC 443 (DeNoyer 1979a). Turner (1969) suggested that the 
inversion of the 1720 MHz line is produced by infrared pumping. Since then, however, 
galactic OH maser emission at 1720 MHz near SNRs has been largely forgotten. In the 
following years, the evidence that the three SNRs above are interacting with nearby clouds 
has grown as to make them three of the best known and studied cases of SNR-molecular 
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cloud interactions (eg. Wootten 1977, 1981, DeNoyer, L.K. 1979b, 1983, White et al. 1987, 
Burton et al. 1988). It is natural then to suggest that this interaction plays a key role in 
the production of maser emission at 1720 MHz (Frail et al. 1996). Several OH 1720 MHz 
surveys have been carried out recently (Frail et al. 1996, Green et al. 1997, Koralesky et 
al. 1998) in order to find new cases of SNR-OH 1720 MHz maser associations. At present, 
among 150 surveyed Galactic SNRs, positive results have been obtained for around twenty 
of them. 

The origin of this maser emission can be easily explained through collisional excitation 
by the passage of a shock front (Frail et al. 1996). Therefore, the presence of OH 1720 MHz 
masers is proposed as a powerful tool to diagnose SNR-molecular cloud interactions. Elitzur 
(1976) showed that the pumping is more efficient if the kinetic temperature is in the range 
25 to 200 K, and the H 2 density is in the range 10 3 to 10 5 cm" 3 . Locket, Gauthier & Elitzur 
(1999) found tighter constrains on the physical conditions needed to produce OH 1720 MHz 
maser emission: the presence of masers imply moderate temperatures (50 - 125 K) and 
densities of the order of 10 5 cm -3 , which can exist only if the shocks are of C-type. Thus, 
the abscence of OH 1720 MHz masers in itself is not enough to rule out an interaction. 
Since the maximum amplification occurs along the edge of an SNR, where the least velocity 
dispersion and hence the largest coherence is maintained, the velocity of an OH 1720 MHz 
maser associated to an SNR can be identified with the systemic velocity of the remnant, 
making it possible to constrain its distance. This property is also suggested by Claussen et 
al. (1997), who point out that the only OH 1720 MHz maser emission observed toward IC 
443 appears at the site where clear evidence for the existence of a transverse shock is found. 

In a recent paper, Frail & Mitchell (1998) surveyed the regions of the masers associated 
with the well studied SNRs W28, W44 and 3C 391, in high excitation CO and H 2 CO 
transitions. The latter SNR is also known to be interacting with molecular gas (Wilner, 
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Reynolds & Moffett 1998, Reach & Rho 1988, 1999). The authors succeeded in finding the 
molecular shocks in all three of them, supporting the idea that the masers originate in post 
shock gas. They found that the masers are distributed near, but not coincident with, peaks 
of CO J=3-2 emission. 

Motivated by the proposed interaction between SNRs and molecular clouds that OH 
1720 MHz masers may reveal, we have selected three SNRs: G349.7+0.2, G16.7+0.1 
and CTB 37A, and surveyed the molecular gas around them in search for the clouds 
hypothetically interacting with the shock fronts. We used the lowest energy transition of the 
CO molecule as a probe, which is ubiquitous and easy to detect, as a first approach to the 
structure of the ISM around these remnants. The present observations allowed us to detect 
several sites of possible interaction between these SNRs and the surrounding molecular 
gas, even at locations where no OH 1720 MHz masers are present (e.g. G16. 7+0.1). These 
observations provide the basis for future work analyzing molecular shocks through higher 
transitions, and for comparing the physics of shocked and unshocked gas. 

2. Observations and data reduction: 

Observations of the CO J=l-0 line (115.271204 GHz) were performed with the NRAO 2 
12 Meter Telescope located on Kitt Peak, Arizona, on 1998 June 12-18 and 27-30. The 
beam size at this frequency is 54". A square field of ~ 38' x 38' was covered toward each 
of the remnants. Observations were made using the On-The-Fly (OTF) technique, allowing 
efficient imaging of extended regions. Three difrerent spectrometers were employed: two 
of them with 128 channels and resolutions of 1 MHz and 500 kHz, respectively, and 

2 The National Radio Astronomy Observatory is a facility of the National Science 
Foundation operated under cooperative agreement by Associated Universities, Inc. 
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a 768-channel correlator with a resolution of 98 kHz. An absolute position-switching 
mode was used, with reference positions given in Table [l|. These reference positions were 
measured to be free of CO emission down to levels of less than about 0.5 K. Four images 
were obtained for G349. 7+0.2, five for G16. 7+0.1 and two for CTB 37A. Since data at the 
12 Meter Telescope are calibrated by default to the T* R scale, which is equivalent to the 
radiation temperature T# for a source much larger than the main diffraction beam of the 
telescope, it is often necessary to apply corrections which refer to sources of different sizes. 
For sources with sizes in the range 20" — 40", the efficiency factor which converts T^ to Tr 
is 0.85 based on historic measurements of the planets Jupiter and Saturn. 

The AIPS package was employed for data processing. Baselines were subtracted in the 
UV plane directly with the task SDLSF, which uses a unique range of line-free channels 
in the whole area. Because of the proximity of all three sources to the Galactic plane, 
the intensity distribution along the observed images made it extremely difficult to find a 
suitable range of line-free channels. We thus decided to give priority to the baseline removal 
in the regions closest to the remnants and OH 1720 MHz masers, with possible detriment 
to the rest. Images were averaged for each source using the task WTSUM, with 1/cr 2 as 
weight, where a is the noise image corresponding to each image. Due to the very high noise 
of one of the original images of G349.7+0.2, only three of them were kept in the average for 
this source. 

In Table [TJ we summarize the observational parameters, including the noise of the final 
images for each SNR. 



EDITOR: PLACE TABLE | HERE. 
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Table 1. Observational Parameters 



Frequency [GHz]: 


115.27 (2.6 mm) 






Spectral resolution [kHz, km s _1 ] 


98 


0.25 


768 


and number of channels: 


500 


2.6 


128 




1000 


1.3 


128 


Angular resolution [arcsec]: 


53.6x53.6 






Surveyed area [arcmin]: 


~ 38 x 38 






Observing mode : 


absolute position switching 








G349. 7+0.2 


CTB37A 


G16.7+0.1 


Central position [RA (1950)] 


17 14 36.0 


17 11 00.0 


18 18 07.4 


[decl. (1950)]: 


-37 22 59 


-38 30 00 


-14 23 18 


Reference position [RA (1950)] 


17 00 36.0 


16 35 00.0 


18 35 07.4 


[decl. (1950)]: 


-37 22 59 


-38 00 00 


-14 23 18 


Number of images: 


3 


2 


5 


Observing time [hours]: 


4.1 


3.6 


8.3 


Central velocity (LSR) [km s^ 1 ]: 


-40 


-65 


+50 


rms [K]: 


0.6 


0.8 


0.3 
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3. Results: 

In what follows, we will devote a subsection to each SNR, including a brief introduction 
to the remnant, a description of our CO results and an estimation of physical parameters for 
the structures most probably related to the remnant. Distances shall be derived from the 
Galactic rotation model by Fich, Blitz & Stark (1989). To calculate the H 2 column density, 
the relationship X = N(H 2 ) /Wco will be used, where Wco is the integrated CO line 
intensity, Wco — I Tr(CO) dv, in K km s^ 1 , and where we have divided our antenna-based 
T* R temperatures by the efficiency factor of 0.85 to place them on the Tr scale appropriate 
for sources with sizes in the range 20" — 40". Several studies suggest a trend of increasing 
X with Galactocentric distance (Digel et al. 1997), thus the value of X may vary from one 
remnant to another. 

Masses will be estimated by integrating N H2 over the solid angle subtended by the CO 
emission feature. Assuming spherical volumes and using a mean molecular weight per H 2 
molecule of 2.72 m H (Allen 1973), the mass and density of a cloud can be expressed as: 

M = 5.65 x 10~ 21 X Wco 9 2 D 2 M Q (1) 

and 

4.1 x IP" 19 X Wco Vf _ 3 f0 , 

n H2 = cm d , (2) 

where D is the distance in kpc, 9 is the angular radius in arcmin of (the cloud assumed 
spherical), computed as half the measured diameter deconvolved by the beam size, and / is 
a correction factor allowing for elongated shapes, and is given by the ratio 9 ma j/9 m i n . 
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3.1. G349.7+0.2: 

This remnant consists of an incomplete clumpy shell, 2' in diameter, with enhancement 
to the south (Figure 1). HI observations (Caswell et al. 1975) suggest that G349. 7+0.2 
is located beyond the tangent point, since absorption features are detected even at 
positive velocities, up to about +6 km s _1 . Clark & Caswell (1976) notice that the strong 
north-south emission ridge, somewhat east of the center, probably represents a non-uniform 
distribution of emission in a shell rather than a filled center, and would be evidence for a 
strong interaction with the ISM at an early stage in the life of this SNR. Frail et al. (1996) 
detected five OH 1720 MHz masers, with the brightest of them aligned along the emission 
ridge, and a few weak ones at the southern continuum knot. 

Our CO observations clearly reveal a concentration inside the remnant, at 
RA=17 h 18 m O a .O, decl=-37°26'36" (J2000). This concentration appears just above the 
southern knot, coincident with the emission ridge. A gaussian fit to the spectrum of this 
cloud yields a central velocity of +16.5 km s _1 and a FWHM velocity of Av = 4 km s _1 . 
Figure 2 shows the cloud integrated over this Av. There is a striking coincidence among 
the CO cloud and the OH 1720 MHz masers, which appear encircling the former. From 
Table 2 in Frail et al. (1996), the velocities of the OH 1720 MHz masers vary from +14.3 
to +16.9 km s -1 , in excellent agreement with the velocity of the CO cloud. 

If we adopt the systemic velocity of the cloud to be +16.5 km s _1 , the distance turns 
out to be ~ 23 kpc, in close agreement with previous results (Frail et al. 1996). The 
diameter of the cloud is estimated to be ~ 50 arcsec, which corresponds to 5 pc. This CO 
cloud is placed at a Galactocentric distance of R = 14 kpc. For R = 12-13 kpc, Digel, Bally 
& Thaddeus (1990) estimate that X = (8 ± 4) x 10 20 cm' 2 (K km s -1 ) -1 , while Mead k 
Kutner (1988) give a value of ~ 4 cm -2 (K km s -1 ) -1 . Therefore, we use X = 8 x 10 20 
cm" 2 (K km s -1 ) -1 keeping in mind that the mass and density may be over estimated by 
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a factor of 2. With these values, the mass and H 2 number density yield 1.2xl0 4 M Q and 
1100 cm -3 , respectively. 

In Table ^| we summarize the results obtained above. The second column quotes the 
approximate central position of the structure in equatorial coordinates referred to J2000. 
The third column lists the radius r after deconvolution with the beam size. The fourth and 
fifth columns contain the systemic velocity v sys and the FWHM velocity Av, both in km 
s _1 . The computed mass and H 2 density are given in the last two columns. 

EDITOR: PLACE TABLE @ HERE. 
3.2. CTB 37A: 

This source (Fig. 3) is actually two SNRs overlapped in projection (Kassim, Baum & 
Weiler 1991). One of them, G348. 5+0.1 (to the west in Fig. 3), is a partial shell with a 
faint extension out the open end of the shell ("breakout" morphology), while the other one, 
G348. 5-0.0 (to the east in Fig. 3), taken to be a jet of the former by Milne et al. (1979), is 
a partial shell. By comparing HI spectra of CTB 37A with G349.7+0.2, which is about one 
degree apart, Caswell et al. (1975) conclude that CTB 37A is beyond the tangential point 
(at -110 km s _1 ) but not farther away than a feature seen at -65 km s _1 . While most of the 
OH 1720 MHz masers identified by Frail et al. (1996) appear projected over G348. 5+0.1 
between -63.5 and -66.3 km s _1 , two of them fall at -21.4 and -23.3 km s _1 respectively. 
These two latter masers appear in the region where the two remnants overlap, suggesting 
that they are associated with G348. 5-0.0. Therefore, the OH 1720 MHz masers allow an 
assignement of the systemic velocities of the two remnants that constitute CTB 37A. 

An inspection of the CO emission toward CTB 37A throughout the whole observed 
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velocity range, reveals a very complex distribution. In searching for a positional coincidence 
among the OH 1720 MHz masers and the CO emission, the best association could be found 
from -22.7 to -25.3 km s _1 (Fig. 4) and from -60.4 to -68.3 km s" 1 (Fig. 5). In what 
follows we will describe the identified structures likely to be interacting with the CTB 37A 
complex. 

At about -24 km s" 1 , a weak CO concentration is observed to the west of G348. 5-0.0, 
with a few brighter clumps embedded. This structure, centered approximately at 
RA=17 ft 14 m 39 s , decl=-38°30' (J2000), shall be called hereinafter "the western cloud". 
The two OH 1720 MHz masers at -21.4 and -23.3 km s -1 lie on the edge of the cloud. 
A Gaussian fit to the average profile towards the western cloud yields a central velocity 
of -23.1 km s" 1 with a FWHM velocity of 5.3 km s" 1 . At this systemic velocity, the 
kinematical distance could be 3.1 or 13.5 kpc. If we assume that this cloud is related to the 
CTB 37A complex, the HI absorption study of Caswell et al. (1975) solves the ambiguity 
in favor of the larger distance. Regardless of the ambiguity, the Galactocentric distance is 5 
kpc. 

A value of X = 2 x 10 20 cm" 2 (K km s -1 )" 1 will be used in deriving mass and mean 
H2 density, based on 7— ray emission studies for the inner Galaxy (Lebrun et al. 1983). 
Assuming an ellipsoidal geometry with the major and minor axes equal to 340" and 135" 
respectively, the mass is estimated to be ~ 15900 M , and the average H 2 density, 150 
cm -3 . 

Other structures probably associated with the OH 1720 MHz masers are found beyond 
-60 km s -1 . In the integrated image shown in Fig. 5, there is an extended structure covering 
the north-west half of G348.5+0.1. Several components can be recognized inside this large 
structure. Because of the proximity with the OH 1720 MHz masers, we will focus on two of 
them: one centered near RA=17 /l 14 m 20 s .5, decl=-38°33'30" (J2000) with three OH 1720 
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MHz masers on its northern edge, which shall be called "the central cloud" ; and another 
one located to the north, coincident with the continuum shell from RA=17 h 14 m 50 s .5, 
decl=-38°27'10" (J2000) to RA=17 /l 14 m 17 s .5, decl=-38°30'40" (J2000), which shall be 
called "the northern cloud" . Another small, weak cloud which appears detached from the 
bulk of the CO emission will also be studied. This cloud, hereafter called "the southern 
cloud", is centered approximately at RA=17 /l 14 m 28 s .5, decl=-38°41'50" (J2000) and 
contains an OH 1720 MHz maser on its northern edge. 

Since the three clouds are at -65 ±2 km s _1 , we adopt a Galactocentric distance of 
3.5 kpc for them, and a kinematical distance of 11.3 kpc, where again we have discarded 
the closer distance of 5.3 kpc based on the HI absorption results. The central and FWHM 
velocities of each cloud, as well as the derived masses and densities, are listed in Table |3|. 
In this table, we also include the western cloud at -23.1 km s _1 . Units are as in Table 
0. In the third column, the listed radii are derived assuming spherical geometries and, 
where applicable, the correction factor / is included in parenthesis. The northern cloud is 
assumed to be an ellipsoide enclosed by the outermost white contour in Fig. 5. In fitting 
a Gaussian to the profile towards the central cloud, only the highest peak was taken into 
account, since the profile looks highly asymmetric and consisting of several components. 
For all the clouds, X was assumed to be equal to 2 x 10 20 cm" 2 (K km s -1 ) -1 . 

EDITOR: PLACE TABLE | HERE. 
3.3. G16. 7+0.1: 

G16. 7+0.1 (Fig. 6), which belongs to the composite type of SNRs, consists of an 
irregular shell, brighter to the south, and a central jet-like feature probably powered by an 



Table 2. Observed and derived parameters for the CO structure associated with 

G349. 7+0.2 



central coordinates r v sys Av M nn 2 

RA, decl. (J2000) (arcsec ) (km s" 1 ) (km s" 1 ) (M ) ( cm 

17 18 0.0,-37 26 36 24 +16.5 4.0 1.2xl0 4 1100 



Table 3. Observed and derived parameters for CO structures associated with CTB 37A 



cloud 


central coordinates 


r 


Vsys 


Av 


M 


71h 2 




RA, decl. (J2000) 


(arcsec) 


(km s _1 ) 


(km s _1 ) 


(M ) 


(cm' 3 ) 


Western 


17 14 39, -38 30 


107 (2.5) 


-23.1 


5.3 


1.5xl0 4 


150 


Central 


17 14 20.5, -38 33 30 


60 (2.1) 


-64.0 


6.9 


7.2xl0 3 


520 


Northern 


17 14 30, -38 27 


150 (3.13) 


-65.0 


12.0 


5.8xl0 4 


660 


Southern 


17 14 28.5, -38 41 50 


37 


-66.6 


4.0 


1.3xl0 3 


280 
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undetected pulsar (Frail & Moffett 1993). Green et al. (1997) detected a single OH 1720 
MHz maser near the bright southern edge of the shell, at +20 km s _1 . Unfortunately, there 
are no line absortion studies to constrain the systemic velocity of the remnant and compare 
it to the velocity of the OH 1720 MHz maser. 

Like CTB 37A, G16. 7+0.1 appears to be located in a very complex region. The CO 
emission is extended over a large velocity range. The only feature found at the position 
of the OH 1720 MHz maser is a small cloud at ~ +25 km s^ 1 , centered approximately 
at RA=18 /l 20 m 56 s .2, decl=-14 21'55" (J2000). This feature will be called "the southern 
cloud". The 5 km s _1 gap between the velocities of the OH 1720 MHz maser and the 
southern cloud, may be related to the anomalous velocity dispersion of the maser, larger 
than typical values (Green et al. 1997). The southern cloud is seen from +25.1 to +25.9 
km s _1 . Figure 7 shows the CO emission integrated over this velocity range. 

There are two more structures that seem to be associated with the remnant. One 
of them is a small concentration inside the shell, centered at about RA=18 h 20 m 59 s .6, 
decl=— 14°19'22" (J2000), located to the east of the central continuum feature. In what 
follows, we will call this concentration "the central cloud". The other structure is a 
bright concentration to the northwest, hereinafter "the north-western cloud", centered 
approximately at RA=18 h 20 m 46 s .6, decl=-14°18'10" (J2000). The close agreement among 
the contours of the CO and continuum emission, strongly suggests that the north-western 
cloud is causing the flattening of the shock front observed in this direction. 

The parameters computed for these three clouds are listed in Table f|. To estimate the 
distance, we first notice that the peak velocities of the clouds are within 0.2 km s _1 of +25.6 
km s _1 . At +25.6 km s _1 , the Galactic rotation model predicts two kinematical distances: 
2.6 and 13.7 kpc, and a Galactocentric distance of 3.5 kpc. As stated above, there are no 
absorption studies to solve the kinematical distance ambiguity. Therefore, we attempted to 
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find which of the two distances is most likely applicable using a Ti-D relationship, in spite of 
the large intrinsic errors of this method (typical uncertainties can be more than a factor of 
2). Based on the surface density given by Helfand et al. (1989), the H-D curve estimated by 
Case & Bhattacharya (1998) produces a distance of 13.7 kpc, in excellent agreement with 
the farther distance given by the Galactic rotation model. Though this perfect agreement 
may a coincidence, the H-D distance implies that the lower kinematical distance should be 
discarded. The ratio X is assumed to be equal to 2 x 10 20 cm -2 (K km s -1 ) -1 . The units 
are as in Table ^. The third column is as described in Table |3|. 

EDITOR: PLACE TABLE | HERE. 



4. Discussion: 

4.1. Correlation between CO clouds and OH 1720 MHz masers 

To compare quantitatively the agreement between the velocities of the OH 1720 MHz 
masers and the velocities of the CO peak flux densities, CO profiles toward all individual 
masers were taken, and Gaussian fits to the identified peaks were calculated. The average 
difference in velocities is |AV| = 0.6 km s _1 , with a disperssion of 1.7 km s -1 . The highest 
departure, above 5 km s" 1 , corresponds to G16. 7+0.1. Disregarding this latter value, |AV| 
turns out to be 0.2 ± 1.2 km s -1 . 

In addition to the velocity agreement, the positional agreement between CO clouds 
and OH 1720 MHz masers is quite striking in many cases. The best match is found in 
G349.7+0.2, where all of the masers are aligned around the only CO feature detected 
towards the remnant. No correlation is found between maser peak flux densities and 
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proximity to the CO emission peak. Brogan et al. (2000) were able to measure magnetic 
fields of ~ 300/iG using the Zeeman effect towards the brightest OH masers in this SNR. 

Another case in which a very good agreement is found is G16. 7+0.1. The OH 1720 
MHz maser is not only superimposed on a CO clump (the southern cloud) but also lies at 
the contact region between the southern cloud and a continuum maximum. However, the 
velocities are not coincident. According to Green et al. (1997), the maser is at ~ +20 km 
s _1 . The CO images around this velocity show no emission near the position of the OH 
1720 MHz maser. In the next paragraph, some suitable explanations for this discrepancy 
are proposed. 

Claussen et al. (1997) found that some of the OH 1720 MHz masers around the SNR 
W28 do not share the systemic velocity of the remnant, but are blueshifted by ~ 10 km s _1 , 
as is the molecular gas associated with them. The CO profile towards the southern cloud 
in G16. 7+0.1, shows a weak wing centered at +19.3 km s _1 , with a 6.5 km s" 1 FWHM. 
This wing, which is in the limit of detect ability, would imply that part of the molecular gas 
of the southern cloud was accelerated to ~ 6.5 km s" 1 by the shock wave and, like in the 
case of W28, the OH 1720 MHz maser indicate the velocity of the accelerated gas. It is 
also possible that the maser is in fact blended features, as Green et al. (1997) suggested. 
In such case, it would be important to compare the velocities of splitted components to the 
systemic velocity of the southern cloud to analyze the agreement among them. 

In CTB 37A, the correlation is poor compared to the other two SNRs. Although the 
masers are in most cases located inside CO structures, the proximity to emission peaks 
does not always apply. Moreover, the strength of the masers is unrelated to the local CO 
emission intensity. There are two masers, located at RA=17 /l 13 m 54 s .87, decl=— 38°33'25 // .5 
(J2000) and RA=17 fe 13 m 58 s .81, decl=-38°34'26".2 (J2000), with velocities -65.1 and -65.2 
km s™ 1 respectively, which lie outside the edge of the CO emission region. Brogan et 
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al. (2000) found that, unlike W28 and W44 where magnetic fields are uniform both in 
magnitude and direction, in CTB 37A the magnetic field strength B 9 changes direction 
along the line of sight over length scales of only ~ 3 pc, revealing a complicated magnetic 
field morphology. 

The CO features found with the present data to be associated with the OH 1720 
MHz masers, have densities not higher than ~ 10 3 cm" 3 , which is at least two orders of 
magnitude lower than theoretical predictions. The 12 CO J=l-0 line is highly saturated 
toward the dense center of molecular clouds, and thus does not vary linearly with the 
column density (Oka et al. 1998). However, low-density cloud envelopes, which are less 
opaque, are sensitively traced by this line. Since most of the material in a cloud resides 
in the envelopes, this line is adequate to estimate a reliable mass for the cloud. To detect 
high-density clumps or cores inside clouds, with densities in the range 10 4 to 10 6 cm -3 , 
other species like CS, C 18 or H 2 CO (Mangum & Wootten 1993, Bronfman, Nyman & May 
1996) or higher CO transitions (eg. Frail & Mitchell 1998) need to be observed. Frail & 
Mitchell (1998) estimated the densities of post-shock gas associated with OH 1720 MHz 
masers to be n H , 2 ~ 10 4 — 10 6 cm -3 , and infer pre-shock densities of the order 10 2 to 10 3 
cm" 3 . In summary, our data reveal the ISM structures with which SNRs are interacting, 
but do not probe the shocked molecular gas in which the OH 1720 MHz maser emission is 
produced, which should be observed in other lines. 

4.2. Interaction between SNRs and molecular clouds 

Large molecular complexes are proposed to be the birthplace of the young, massive 
progenitor stars of Type II supernovae. Half of the galactic SNRs are supposed to be 
of this type, and this would be the reason why Huang & Thaddeus (1986) found a very 
good correlation between large molecular cloud complexes and half of the SNRs surveyed 
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in CO. It is likely that the observed molecular complexes be part of the natal-clouds of 
the pre-supernova stars. However, although it is statistically improbable that such a good 
correlation arises from chance superpositions, there is very little direct evidence that SNRs 
are actually in physical contact with the molecular clouds observed. 

Shocked-excited 1720 MHz OH masers toward SNRs constitute a powerful tool to 
ensure on safe grounds that the molecular structures at the same velocity of the masers 
are physically related to the remnants. The use of OH 1720 MHz masers as signposts of 
SNR-molecular cloud interactions allowed us to detect the molecular gas into which SN 
shock fronts are being driven. In what follows, we analize the interaction of each SNR with 
its surrounding molecular gas. 

4.2.1. G349.7+0.2 

In this case, the evidence for SNR-molecular cloud interaction is compelling, based 
on the positional and velocity coincidence between the CO feature and the OH 1720 
MHz masers. This conclusion supports the statement of Clark & Caswell (1976) that 
the expansion has been influenced by a dense environment in early stages. To obtain an 
additional proof of interaction, we have searched for IRAS point sources that fulfill the 
following conditions: (1) S 100 > 20 Jy, (2) 1.2 < S wo /S m < 6.0, and (3) S 60 /S 25 > 1, 
where S\ denotes the (uncorrected) IR-flux density in the wavelength A /im. These selection 
criteria are indicative of proto-stars (Junkes, Fiirst & Reich 1992) or dust heated in SNR 
shocks (Arendt 1989). We have found that the point source IRAS 17146-3723, which agrees 
with these conditions, is positionally coincident with the CO cloud. It is unlikely that this 
IRAS source is a protostar triggered by the shock wave of G349.7+0.2 because of the time 
scales involved. Formation of protostars from condensed molecular matter needs a period of 
the order of 10 6 years, enough for the matter to cool, which is too long compared to the age 
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of a SNR. We therefore suggest that IRAS 17146-3723 is produced by shock-heated dust 
and, together with the OH 1720 MHz maser emission, confirms the interaction between the 
SNR and the molecular cloud detected here. A thorough analysis of the more extended 
region surveyed here, is given in Reynoso & Mangum (2001). 

4.2.2. CTB 37 A 

Our observations show two distinct clouds towards CTB 37A, both coincident spatially 
and in velocity with the two groups of OH 1720 MHz masers. Thus, the present CO 
data provide additional support to the theory of the two overlapped SNRs. However, the 
complexity of the CO distribution in this direction makes it difficult to identify molecular 
clouds interacting with the SNR shock fronts. One of the remnants, G348.5+0.1, has 
a characteristic "breakout" morphology, which would imply a strong density gradient. 
Nevertheless, there is no correspondence between the CO distribution and such radio 
morphology. This is in contrast with the case of 3C 391, where the 12 CO distribution 
shows a molecular cloud exactly parallel to the bright radio emission ridge (Wilner et al. 
1998). As in the case of G349. 7+0.2, an additional evidence for interaction is given by the 
IRAS point source 17111-3824, which has an IR spectrum typical of shocked heated dust 
according to the criteria of Junkes et al. (1992), and lies close to the OH maser at the 
northern cloud. This IRAS source is shown as a white star in Fig. 5. 

Apart from the presence of the OH 1720 MHz masers and IRAS 17111-3824, there are 
no further indications of shocked molecular gas such as line broadenings or wings. However, 
at ~ —88 km s _1 there appears a bright, elongated CO structure positionally following the 
masers closest to the central cloud. In Fig. 8, two profiles are shown, the upper towards 
the three masers at the central cloud (white crosses in Fig. 5), and the lower towards an 
arbitrary direction displaced from the remnant, at the northwest of the field. In the upper 
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profile, three peaks are cleary distinguished: one at ~ —65 km s" 1 , which corresponds to 
the gas in which the OH 1720 MHz masers originate, and the others at approximately 
-88 and -107 km s _1 . The peak at ~ —107 km s^ 1 is also seen in the profile outside the 
SNR, and is commonly regarded as an expanding arm about 3 kpc from the Galactic center 
(Caswell et al. 1975). The other peak, instead, does not appear in the lower profile, so it is 
likely to be related with the central cloud. At -88 km s _1 , the only CO emission observed 
extends over the area enclosed by the radio continuum emission. Therefore, we propose 
that the CO emission at this velocity corresponds to a fraction of molecular gas accelerated 
by the shock front towards the rear side of the remnant at a velocity of about 25 km s _1 . 
This accelerated component of the gas would be additional proof for an SNR-molecular 
cloud interaction. 

Chevalier (1999) proposes that in the initial interaction, the (radiative) SNR shell 
drives a slab into the cloud. We suggest that in this case, the slab can be seen as the -88 km 
s _1 component. The post-shock pressure is p c iv 2 d ~8x 10~ 9 dyn cm -2 , where the subindex 
cl indicates that the parameters are estimated at the cloud. Assuming constant density for 
the radiative shell, and allowing for magnetic fields, Chevalier's model states that the ratio 
between the ram pressure in the slab and that at the front of the radiative shell, PdV d / Po^si 
is limited to the range ~ 10 — 100. In this equation, p is the intercloud density and v rs is 
the velocity of the radiative shell before entering the cloud. Thus, assuming n = 15 cm~ 3 
as a typical value for the intercloud density (e.g. Chevalier 1999), the shock velocity at the 
cloud is ~ 15 — 30 km s" 1 . At velocities lower than ~ 50 km s~ x , shocks are of C-type 
(Neufeld & Dalgarno 1989), thus this estimate is consistent with the theoretical prediction 
by Lockett et al. (1999) that the conditions under which 1720 MHz OH maser emission is 
produced can exist only if shocks are of C-type. An upper limit of v rs = 60 km s _1 can be 
reached if an intercloud density as low as n = 1 cm~ 3 (e.g. McKee & Cowie 1975, Dubner 
et al. 1999) is assumed. 



4-2.3. G16. 7+0.1 



Like G349.7+0.2, this is another case in which the evidence for SNR-molecular cloud 
interaction is clear. What is remarkable about this case is that one single maser was found 
to be associated with a region with several structures of molecular gas and more than one 
site of possible interaction between this gas and the SNR shock front. The same occurs 
with IC 443 which, in spite of being one of the remnants more clearly interacting with dense 
molecular gas, has only one maser coincident with a site where the shock is transverse 
(Claussen et al. 1997). 

Based on high-resolution observations, Frail & Moffett (1993) describe the central 
jetlike radio structure as consisting of a center of symmetry at RA=18 fe 20 m 57 s .48, 
decl=— 14°19 / 57".6 (J2000) from which emission extended in two directions opens up into 
more extended lobes. In Fig. 7, it can readily be seen that the northern lobe seems to be 
inclined towards the central cloud, tending to enclose the inner CO emission peak (see also 
Fig. 6). Also, the southern lobe appears to the east of the southern cloud, with the OH 
1720 MHz maser lying at the contact region between them. Both lobes show enhancements 
in the radio continuum emission at the locations of the molecular clouds, which is an 
indication of interaction. 

In virtue of the morphological agreement, the north-western cloud appears as an 
excellent candidate for interaction with this SNR, although no line broadenings are detected. 
The interaction between G16. 7+0.1 and the northern cloud places this remnant in the 
same class as IC 443 and 3C 391 in the sense that, although the morphologically most 
convincing site of SNR-molecular gas interaction is suggested by a flat radio continuum 
rim exactly parallel to the boundary of a molecular cloud, there appears only one maser 
(actually two in the case of 3C 391), and not at that site but near the opposite side of the 
remnant. Moreover, the morphology involved in such cases makes it difficult not to conceive 
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a strongly transverse shock. Thus, it appears that not all transversely-shocked dense gas 
regions lead to the conditions which give rise to OH 1720 MHz maser emission. Possible 
explanations for the lack of OH 1720 MHz masers in these cases could be related, for 
instance, to dissociation of OH by high velocity shocks. The trend observed in these three 
remnants is undoubtfully a topic that deserves further investigation, and should certainly be 
taken into account in future theoretical works. 

5. Conclusions: 

The surroundings of the SNRs G349.7+0.2, CTB 37A and G16.7+0.1 have been 
investigated in the CO J=l-0 transition. These remnants have been selected since all 
of them have associated OH 1720 MHz maser emission. Our observations allowed us to 
determine the density distribution of the ISM into which these remnants are expanding, 
revealing a rich environment in all three cases. In spite of the confusion produced by 
the ubiquity of the CO J=l-0 emission, we have been able to identify molecular clouds 
interacting with SNRs based mainly on the positional and kinematical coincidence with 
shock-excited OH 1720 MHz masers, but also on morphological arguments as well as 
radio continuum enhancements, IR emission by shock- heated dust, and peculiar spectral 
components. Unfortunately, we have not detected spectral line broadenings, which are 
unmistakable signatures of interaction. However, there are other cases of SNRs interacting 
with molecular clouds, in which no line broadenings have been observed in the CO J=l-0 
transition (Wilner et al. 1998, Dubner et al. 1999). Reasons for this non-detection vary 
from rapid dissociation of the molecules, masking of weak CO J=l-0 features by background 
cloud emission, and beam dilution of highly localized blast-wave interactions (Wilner et al. 
1998). 

The main results of this paper are summarized below. 
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1) All but two of the OH 1720 MHz masers are associated with CO emission peaks 
that posses an excellent velocity correlation with the former. The difference between both 
velocities is in average (0.6 ± 1.7) km s" 1 in the worst case. 

2) The masers are positionally close to, but not coincident with, CO peaks. The same 
trend has been observed by Frail & Mitchell (1998). 

3) The CO emission associated with the OH 1720 MHz masers allows constraint of the 
systemic velocities of the SNRs. 

4) The H 2 densities inferred (< 10 3 cm" 3 ) are representative of the preshock gas, but 
do not trace high density clumps or shocked, compressed material. 

5) Several cases of SNR-molecular cloud interaction have been detected based on wings, 
shifted spectral components, morphological coincidences and IR emission from shock-heated 
dust. 

6) In G16. 7+0.1, the OH 1720 MHz maser seems to originate in molecular gas 
accelerated to 6.5 km s -1 , and not in the gas located at the systemic velocity of the SNR. 

7) The interstellar gas does not account for the breakout morphology of G348. 5+0.1, 
one of the SNRs of the CTB 37A complex. 

8) In CTB 37A, the shock is computed to be of C-type, in agreement with theoretical 
predictions (Lockett et al. 1999). 

9) G16.7+0.1 is found to belong to a group of SNRs (together with 3C 391 and IC 443) 
with the following common characteristic: part of the rim is obviously flattened by a dense 
molecular cloud. The morphology strongly suggests that the expansion is perpendicular to 
the line of sight at this location. However, no OH maser emission at 1720 MHz is detected. 
The only OH 1720 MHz masers observed in such remnants (only one or two) are far from 
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the flattened part of the rim. This behavior seems to contradict the suggestion that OH 
1720 MHz masers appear where shocks are transverse. 

We are grateful to the NRAO for allocating time in the 12 Meter Telescope for 
this project. We also thank D. Helfand, N. Kassim and C. Salter for providing us their 
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(Argentina) for visiting the 12 Meter Telecope, as well as NRAO support and staff assistence 
during her visit. This research was partially funded through a Cooperative Science Program 
between CONICET and NSF and through CONICET grant 4203/96. 
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Table 4. Observed and derived parameters for CO structures associated with G16. 7+0.1 



cloud 


central coordinates 


r 


Vsys 


Av 


M 


n H2 




RA, decl. (J2000) 


(arcsec) 


(km s -1 ) 


(km s -1 ) 


(M©) 


(cm~ 3 ) 


Southern 


18 20 56.2, -14 21 55 


32 (1.6) 


+25.8 


1.5 


2.1 x 10 3 


260 


Central 


18 20 59.6, -14 19 22 


28 (1.5) 


+25.4 


4.4 


1.8xl0 3 


330 


North-western 


18 20 46.6, -14 18 10 


37 (2) 


+26.6 


2.7 


3.4xl0 3 


280 
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6. Figure captions 

Fig. I. — A radio continuum image of G349.7+0.2 at 1.4 GHz obtained with the VLA, kindly 
provided by C. Salter. The grayscale is indicated in mJy/beam on top of the image, while 
the contours are expressed as a percentage of the peak flux density (0.48 Jy/beam) at 1, 2, 
3, 5, 10, 20, 30, 40, 50, 60, 70, 80, 90, and 99%. For better contrast, white contours are 
used over dark shadows. The white crosses show the positions of the OH 1720 MHz masers 
detected by Frail et al. (1996). The beam, indicated in the bottom left corner, is 19". 4 x 4". 5, 
P.A.=-22°.2. 

Fig. 2. — CO emission integrated between +14.5 and +18.5 km s _1 toward G349.7+0.2, in 
greyscale and contours. Data were taken with the 12 Meter Telescope of the NRAO. The 
greyscale is in K km s -1 , and is shown on top of the image. The white contours represent 
the CO emission between 8 and 28 K km s _1 , in steps of 4 K km s _1 . The crosses show the 
positions of the OH 1720 MHz masers. The beam, 60" x 60", is indicated by the open circle 
in the top left corner. The 3a noise level is 6 K km s~\ A few radio continuum contours in 
grey are included to represent the SNR. 

Fig. 3. — A radio continuum image of CTB 37A at 1.4 GHz, from Kassim et al. (1991), 
obtained with the VLA. Two sources can be identified: G348. 7+0.1 to the west, and G348.5- 
0.0 to the east. The grayscale is indicated in mJy/beam on top of the image, while the 
contours are at 1, 3, 5, 10, 20, 30, 65, and 95 % of the peak flux. Here again, white contours 
are used over dark shadows. The peak flux is 0.43 Jy/beam. Black crosses show the positions 
of the OH 1720 MHz masers between -63 and -67 km s" 1 , while white crosses represent the 
two masers at -25.3 and -22.7 km s -1 (Frail et al. 1996). The beam, indicated in the bottom 
left corner, is 32".9 x 18", P.A .=4° .8. 

Fig. 4. — CO emission integrated between -25.3 and -22.7 km s" 1 toward CTB 37A. Data 
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were taken with the 12 Meter Telescope of the NRAO. The greyscale is in K km s" 1 , and is 
shown on top of the image. The beam, 53". 6 x 53". 6, is plotted in the bottom left corner. A 
few representative contours of the continuum emission from CTB 37A are included as grey 
lines. CO emission contours at 11.5, 17, 22 and 28.5 K km s -1 are plotted as black lines, 
except for the last one which appears in white. The two white plus signs represent the two 
OH 1720 MHz masers at -21.4 and -23.3 km s" 1 (Frail et al. 1996). The 3a noise level is 9 
K km s _1 . 

Fig. 5. — CO emission integrated between -68.3 and -60.4 km s -1 toward CTB 37A. Data 
were taken with the 12 Meter Telescope of the NRAO. The greyscale is in K km s -1 , and 
is shown on top of the image. The beam size is 60" x 60". A few representative contours 
of the continuum emission from CTB 37A are included as thick grey lines. CO emission 
contours at 20, 30, 40, 55, 75, and 95 % of the peak emission are plotted in black or white 
solid lines, according to the darkness of the background. The peak emission is 82.2 K km 
s _1 . The plus signs show the positions of the OH 1720 MHz masers found by Frail et al. 
(1996) between -63 and -67 km s -1 , where white lines have been used over dark background 
for better contrast. The white star represents the IRAS point source 17111-3824 discussed 
in the text (Section 4.2.2). The arrows indicate the three clouds discussed in the text. The 
3c noise level is 7 K km s _1 . 

Fig. 6.— A radio continuum image of G16. 7+0.1 at 4.86 GHz from Helfand et al. (1989), 
obtained with the VLA. The grayscale is indicated in mJy/beam on top of the image, while 
the contours are at 1, 5, 15, and 30 mJy/beam. The plus sign indicates the OH 1720 MHz 
maser detected by Green et al. (1997). The beam, indicated in the bottom left corner, is 
21" x 16", P.A.=11°. 

Fig. 7. — CO emission integrated between -25.1 and -25.9 km s _1 toward G16. 7+0.1. Data 
were taken with the 12 Meter Telescope of the NRAO. The greyscale is in K km s _1 , and 
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is shown on top of the image. The beam size, 53". 6 x 53". 6, is plotted in the bottom 
right corner. A few representative contours of the continuum emission from G16. 7+0.1 are 
included as dark-grey lines. CO emission contours at 3.8, 7.6, 12.7, 16.5, 20.3, and 24 K km 
s _1 are plotted in black and white solid lines, according to the darkness of the background. 
The plus sign in white shows the position of the OH 1720 MHz maser found by Green et al. 
(1997). The arrows indicate the three clouds discussed in the text. The 3<r noise level is 1.5 
K km s _1 . 

Fig. 8. — CO emission profiles toward the OH masers lying on the central cloud in CTB 37A 
(upper profile) and towards an arbitrary direction, to the northwest of the observed field 
and displaced from the continuum emission (lower profile). Brightness temperatures are in 
Kelvins, while velocities are in km s _1 . 
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